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Abstract. High resolution diffraction grating spectra of X-ray emission from massive stars ob- 
tained with Chandra and XMM-Newton have revolutionized our understanding of their powerful, 
radiation-driven winds. Emission line shapes and line ratios provide diagnostics on a number of key 
wind parameters. Modeling of resolved emission line velocity profiles allows us to derive indepen- 
dent constraints on stellar mass-loss rates, leading to downward revisions of a factor of a few from 
previous measurements. Line ratios in He-like ions strongly constrain the spatial distribution of X- 
ray emitting plasma, confirming the expectations of radiation hydrodynamic simulations that X-ray 
emission begins moderately close to the stellar surface and extends throughout the wind. Some out- 
standing questions remain, including the possibility of large optical depths in resonance lines, which 
is hinted at by differences in line shapes of resonance and intercombination lines from the same ion. 
Resonance scattering leads to nontrivial radiative transfer effects, and modeling it allows us to place 
constraints on shock size, density, and velocity structure. 
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INTRODUCTION 

Massive stars are of great importance in the cosmic cycle of baryonic matter: ejecta from 
core collapse supemovae are one of the main sources of enrichment of heavy elements 
in the universe; radiation from massive stars, which are the first stars to appear in a star 
formation event, has a profound effect on the surrounding interstellar medium (ISM) and 
thus on the entire star formation process; and their powerful winds further influence the 
ISM through input of kinetic energy, and in later evolutionary stages, enriched material 
from the former core exposed by mass-loss in earlier stages. The winds of massive stars 
are also important for their effect on the evolution of the star itself. The fractional mass 
loss over the core hydrogen burning lifetime of the star can be significant, and in later 
evolutionary stages it is even greater. 

The winds of massive stars are driven by radiation pressure in spectral lines, which is 
a result of their enormous luminosity (10 5 — 10 6 L Q ). The effectiveness of line driving is 
multiplied in the highly supersonic flow as saturated absorption lines are Doppler shifted 
out of their own shadow in the acceleration zone of the wind, allowing for an increase in 
momentum transfer in optically thick lines by a factor of Vco/v t / 7 . The winds of massive 
stars are thus unique laboratories for studies of radiation hydrodynamics. 

Radiative line driving is unstable to perturbations. Although the instability is damped 
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near the base of the wind by the scattered radiation field, hydrodynamic simulations 
show that the instability leads to extreme inhomogeneities in the wind (clumping), as 
well as energetic shocks. This results in significant X-ray emission, which has great 
diagnostic utility in studies of massive star winds. 

The single most important measurable parameter of massive star winds is the mass- 
loss rate, since this sets the scale of evolutionary and ISM feedback effects. Although 
there are many observational diagnostics of mass loss (e.g. Ha recombination lines, 
radio free-free emission, UV P Cygni absorption profiles [1, 2]), all of them are depen- 
dent on modeling assumptions, and thus potentially subject to systematic errors. X-ray 
emission line profiles observed using high-resolution diffraction grating spectrometers 
provide an alternative mass-loss rate diagnostic that does not share the same systematic 
errors as other techniques. X-ray line profiles are insensitive to density inhomogeneity, 
non-monotonic velocity fields, and the ionization of the bulk of the wind. 


X-RAY LINE PROFILE MODELING 

Owocki and Cohen [3] have described a formalism for calculating massive star X-ray 
emission line profiles. The wind is treated as a two-component fluid, with a small frac- 
tion of the wind heated to X-ray emitting temperatures and a dominant cool compo- 
nent which absorbs X-rays. X-rays emitted from the near hemisphere of the wind are 
blueshifted and suffer less absorption, while X-rays emitted from the back hemisphere 
are redshifted and more absorbed. The degree of absorption is characterized by the pa- 
rameter t* = kM / 4jri?* v-x , where x(A) is the microscopic opacity of the wind to X-rays, 
R * is the stellar radius, Voo is the wind terminal velocity, and M is the stellar mass-loss 
rate. This can also be written as r* = »c2*, with 2* = M/AnR^v x the characteristic wind 
mass column density. Example model profiles are shown in the left panel of Figure 1 . 
t* is varied from 0 - 100, with increasing characteristic optical depth corresponding to 
more blueward skewed profiles. This variation in profile shape with wind optical depth 
is the basis for the use of X-ray line profiles as a mass-loss rate diagnostic. The right 
panels of Figure 1 show an example of a model fit to the Fe XVII 15.014 A line of the 
prototypical O star 'Q Pup. The statistical uncertainty on the derived value of t* for this 
line is of order 10%, showing that even data of moderate statistical quality can provide a 
powerful mass-loss rate diagnostic, provided that systematic effects are well understood. 


HE-LIKE TRIPLET LINE RATIOS 

Fine ratios of He-like ions have been used extensively in laboratory and astrophysical 
plasmas as density and UV field diagnostics [4, 5]. The lowest lying excited state, 
Isis : \S) , is metastable, so that collisional or photoexcitation to the 1 si p 3 P levels can 
compete with decay to ground. This alters the observed ratio of the intercombination 
(Is 2 1 So — Islp 3 P; x,y, or alternately i) to forbidden (Is 2 1 Sq — 1 sis 3 S \ ; z, or alternately 
/) line strengths. Because the decay of the ls2s \Sj state to ground is an M2 transition, 
the rate scales with Z 10 , so that by observing He-like ions over a range of a factor of 
two in atomic number (nitrogen to silicon), we probe three orders of magnitude in decay 
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FIGURE 1 . Left panel: comparison of line profiles with different characteristic optical depths r* . Right 
panel: Chandra MEG (top) and HEG (bottom) data for Fe XVII 15.014 A fit with a line profile model 
(red curve) [7] . 


rate. 

Plasma densities in O star winds are too low to affect ffi ratios, but the photospheric 
UV field has a strong effect. The UV field strength, and thus the f/i ratio depends on 
distance from the star. A calculation of this effect is shown for different He-like ions in 
the wind of the early O star £ Pup in the left panel of Figure 2. It is not possible to invert 
the observed f/i ratio to obtain a localized radius of plasma formation, simply because 
the X-ray emitting plasma is spread out over a large range in radii. However, the correct 
radial dependence of the line ratio can easily be incorporated into a line profile model, 
yielding a model that predicts both line shape and line ratios with no additional free 
parameters [6]. An example of such a model fit is given in the right panel of Figure 2. 
Such models have been applied to X-ray spectra of O stars, resulting in good agreement. 


MASS LOSS RATE DIAGNOSTICS 

By measuring the characteristic optical depth r* for every reasonably strong emission 
line in the X-ray spectrum of an O star, it is possible to make a precise measurement 
of the mass-loss rate. We have performed such a measurement for £ Pup, shown in 
Figure 3 [7]. We plot the measured values of t* as a function of wavelength, together 
with two model optical depths with different mass-loss rates. The model optical depths 
are calculated using the equation r* = fc(A)2*, taking the mass-loss rate as the only 
changeable parameter. The dashed line shows the model corresponding to the upper- 
limit mass-loss rate inferred from radio free-free emission, while the solid line shows 
the model corresponding to the best fit to the X-ray data. The dominant uncertainty in 
our mass-loss rate determination is from uncertainty in the microscopic opacity model, 
which results from uncertainties in the elemental abundances and the ionization of 
helium. 
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(Pup 



FIGURE 2. Left panel: dependence of /// ratio in He-like triplets on distance from the photosphere of 
£ Pup [6]. Right panel: Emission lines of Mg XI in the spectrum of c, Pup, fit with a line profile model 
including the radial dependence of/// ratio [6]. Rest wavelengths are indicated with vertical dashed lines. 
From left to right, the lines are resonance (vv), intercombination ( v) , and forbidden (z). 



FIGURE 3. Measured characteristic optical depth r* for strong emission lines in the spectrum of 'C, 
Pup plotted as a function of wavelength [7], The shape of the model curves reflects the microscopic 
opacity of the wind. The dashed and solid lines correspond to mass-loss rates of 8.3 x lO _6 M 0 yr _1 and 
3.5 x lO _6 M 0 yr _1 , respectively. 


RESONANCE SCATTERING 

The resonance and intercombination lines from a He-like triplet are expected to have the 
same profile shape, since they both come from the same ion. Surprisingly, these two lines 
have been found to have different profile shapes for He-like triplets in some O star X-ray 
spectra. In the left panel of Figure 4, we show a model fit to the N VI triplet of 'C, Pup 
[8]. The resonance line is more symmetric than the model, while the intercombination 
line is more asymmetric. This cannot be explained by the radial dependence of the f /i 
ratio, which is included in the model. 

If strong resonance lines are optically thick, the radiative transfer depends on the 
velocity structure along the line of sight. This problem has been well studied in the 
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FIGURE 4. Fits to N VI He-like triplet of £ Pup. A model with no resonance scattering (left panel) 
cannot fit the resonance and intercombination lines simultaneously, while a model including the effects of 
resonance scattering (right panel) fits well. 


context of UV line profile modeling in O star winds. Radial photon escape is due to 
radial velocity gradients ( dv/dr ), while lateral photon escape is due to the spherical 
divergence of the wind (v/r), which results in lateral line-of-sight velocity gradients. In 
typical O star winds, lateral photon escape is favored for r>2R*. Since lateral photon 
escape results in low observed Doppler shifts, strong resonance scattering in X-ray line 
formation results in more symmetric lines. This symmetrization effect only applies to 
resonance lines, which have large oscillator strengths, and not to intercombination or 
forbidden lines. 

This effect can be seen in the models shown in Figure 5 . The strength of resonance 
scattering is described by the characteristic Sobolev optical depth, To,*. Larger values 
of To,* result in more symmetric profiles. This model has been fit to the N VI triplet of 
£ Pup with To,* as a free parameter, as shown in the right panel of Figure 4. The fit is 
greatly improved, indicating that resonance scattering is important in the formation of 
this profile. Evidence for this effect has been found in X-ray spectra of other O stars as 
well. However, detailed interpretation of resonance scattering effects are handicapped 
by the fact that we can only guess at the velocity structure in the X-ray emitting regions. 


FE XVII LINE RATIOS 

The X-ray spectrum of £ Pup shows an anomaly which has not been observed in any 
other astrophy sical X-ray spectrum to date . The ratio of 35 — 2p lines ( 1 7 .05 1 A + 1 7 .096 
A) / (16.780 A) is typically found to be ~ 2.4 for a wide range of astrophysical sources, 
including coronal stars with a wide range of temperatures, and including all other O 
stars. However, in £ Pup, this line ratio is ~ 1.4. There is no obvious mechanism that 
would produce this effect. Photoexcitation of the metastable lowest lying excited state 
(7 = 2) has been considered in a study of Fe XVII line ratios in the magnetic cataclysmic 
variable system EX Hya [9] ; this cannot be effective in this case, because photoexcitation 
is most effective at connecting the upper levels of the 17.051 A and 17.096 A lines, 
respectively (via a manifold of 3 p states). Resonance scattering is also probably not 
important for any of these lines, since none of them has a large oscillator strength. 
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r, = 1 0, q = 0, R 0 = 1.5 
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FIGURE 5. Model line profiles (flux as a function of scaled velocity x) showing the effect of resonance 
scattering. Larger values of the characteristic Sobolev optical depth (to,*) produce more symmetric 
profiles. 


Resonance scattering also has little or no effect on line ratios. 

One possible explanation for the observed line ratio is a near coincidence between 
either the 17.051 A or 17.096 A line and an inner-shell 2p-3d absorption line of a 
low charge state of iron with a large ion fraction in the bulk of the wind, such as 
Fe VI. Such inner-shell absorption lines have been observed in AGN warm absorber 
absorption spectra [10] , and have been the subject of extensive theoretical and laboratory 
investigation [11, 12]; however, lower charge states have received less experimental 
attention and pose greater theoretical difficulties due to the larger number of electrons. 
If such a coincidence exists, then line emission from the X-ray emitting part of the 
wind may be efficiently destroyed in the cool bulk of the wind. This hypothesis may 
also explain why this effect is only observed in £ Pup: if the charge balance of iron is 
significantly different for £ Pup than for the other O stars that have been observed by 
XMM and Chandra, as might be expected from its higher effective temperature, then 
this absorption may only be important in £ Pup. 
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